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Abstract 

A brief review of the experimental situation con- 
cerning the electrically-charged charmoniumlike me- 
son candidates, Z~, is presented. 

Introduction 

The Belle Collaboration reported peaks in the 
iP'tt~ I Q and Xci^~ invariant mass distributions |2l in 
B — s> ip'ir~K [3] and B ->■ Xci7r~.K" EL respectively. 
If these peaks are meson resonances, they would have 
a minimal quark substructure of ccdu and be unmis- 
takeably exotic. However, even though the Belle sig- 
nals have more than 5a statistical significance, the ex- 
perimental situation remains uncertain in that none of 
these peaks have yet been confirmed by other exper- 
iments. An analysis by the BaBar Collaboration of 
B — > rp'ir~K neither confirms nor contradicts the 
Belle claim for the Z(4430)~ ->■ tP'tt' 0. In the 
BaBar analysis, B — > J/ipir~K decays were also 
studied, and no evidence for Z(4430)~ — > J/ip7T~ 
was found. 

In this paper, we review and compare Belle and 
BaBar results on searches for charged charmonium- 
like states. An abbreviated version of this review is 
contained in Ref. (H. 

The Belle observation 

Belle's original Z(4430)~ signal O is the sharp 
peak in the tp'Tr" invariant mass distribution from 
B — >■ iP'tt~K decays shown in Fig.Q] The Belle anal- 
ysis is based on a data sample that is equivalent to an 
integrated luminosity of 605 fb _1 . Figure |2] shows the 
Dalitz plot for B — > ijj'ir ~K candidates, where verti- 
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cal bands for K*(892)^Kir- and K^(U30)-^Ktt~ 
are evident and the Z(4430)~ shows up as a hori- 
zontal band of events between M 2 , _ = 19 and 20 

ip'n 

GeV 2 /c 4 . (In the t/j'-ir^ invariant mass distribution of 
Fig. [U the the K * bands are suppressed by cuts on the 
Kir~ masses ||3l). 

A binned maximum-likelihood fit to the tp'-K" mass 
distribution using a Breit Wigner (BW) resonance 
function for the signal and a polynomial background 
gives a mass of M = 4433 ± 4 ± 2 MeV/c 2 and total 
width of T = 45 1 \l til MeV ' with estimated statis- 
tical significance of more than 6a. Consistent signals 
are seen in various subsets of the data: i.e. for both 
the ip' ->■ l + £~ (£ = e, or fi) and -ip' -4 J/4)Tr + Tr~ 
{J/ip -4 £ + £~) subsamples 0. 

A reflection from Kn dynamics? 

Identifying resonant structures in the tp'ir^ channel 
in three-body B — > ip'ir~K decays has to be done 
with care because of the possibility that dynamics in 
the Ktt~ channel can cause mass structures in the 
ip'ir~ invariant mass distribution that have no relation 
to ip'ir~ dynamics. Energy-momentum conservation 
imposes a tight correlation between the decay angle 
(6 n -) in the Kir~ system [7] and the ip'ir~ invari- 
ant mass, and this results in M 2 ,^ being very nearly 
proportional to cos6 n -. Consequently, interference 
between different partial waves in the Kir~ system 
can produce peaks in Mm^- that are purely "reflec- 
tions" of structures in cos 9 n - and have nothing to do 
with ip'-K~ dynamics. However, in the kinematically 
allowed Ktt~ mass range for B — > ip'ir~K decay, 
only S-, P- and D-waves are significant, and this lim- 
ited set of partial waves can only produce fake ip'ir~ 
mass peaks at a discrete set of mass values. 

In the case of the Z(4430)~, the ip'Ti~ peak mass 
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Figure 1: From Belle [3], the ip'ir~ invariant mass 
distribution for B — > ip'ir'K decays after applying a 
veto on the if* (892) and #1(1430). The open his- 
togram shows the data, and the shaded histogram rep- 
resents the scaled results from the background. The 
solid curves are the fit results. 



corresponds to cos8 w - ~ 0.25, and it is not possi- 
ble to produce a peak near this cos 8 n - value with 
any combination of interfering L =0, 1, and 2 
partial waves without introducing additional, larger 
structures at other cos^- values. This is illus- 
trated in Fig. [3l where the histogram shows the 
distribution of cos 8 n - values for a MC sample of 
B -> Z(U30)-K, Z(U30)-->ip f ir- events where 
the Z(4430)~ mass and width closely correspond to 
Belle's reported values 0. The curves in the fig- 
ure show the results of trying to make a peak at the 
same location with interfering S-, P- and D-partial 
waves in the Kir~ channel. (Here both longitudi- 
nally and transversely polarized tj/'s are considered, 
and no attempt is made to restrict the strength of each 
term to that seen for the 5-, P- and D-wave Kir~ 
components in the data.) These curves show that al- 
though a peak can be made at cos8 n - ~ 0.25, it is 
necessarily accompanied by large enhancements near 
cos8 n - ~ ±1. No such structures are evident in the 
ip'ir~ mass plot of Fig. [T] 
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Figure 2: From Belle [3], the Mjl _ (horizontal) vs. 
(vertical) Dalitz plot distribution for candidate 

B — > ifj'ir~ K events. 

The BaBar search for the Z(4430)" 

BaBar searched HI for the Z(4430)~ in the de- 
cay modes B — > ipir~K, where ip = J/ip or tp' and 
K = K + or Kg. The search was performed using 
the full BaBar data set equivalent to integrated lu- 
minosity of 413 fb _1 (455 xlO 6 BB candidates) col- 
lected at a center-of-mass energy of 10.58 GeV. The 
BaBar search covered two additional decay modes 
(B — » J/t/jTr~K) that were not considered in the 
Belle study. The J/ip modes contain about six times 
(see Table II in Ref. [5]) more data than those with 
the ip' and provide a powerful, high statistics tool for 
understanding the background. The Dalitz plots for 
the four decay modes are shown in Fig. [4] The Ktt~ 
mass distribution shows a clear K*(892) signal in all 
decay modes, and a visible if,* (1430) signal in the 
modes with a J/ip but less significant in the modes 
with tp'. Since the Kn~ system dominates the Dalitz 
plot, the BaBar Collaboration studied the Kir~ mass 
distribution and the angular dependence. Similar con- 
tributions from the different wave intensities are ob- 
tained in the modes with Kg and K + (see Table III in 
Ref. [5]) and, therefore, the modes with K s and K + 
are combined. The Kir~ mass distributions are fitted 
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Figure 3: From Ref. |8], the histogram shows the 
cos 6 n - distribution for a MC -generated ip'n~ reso- 
nance with M = 4.43 GeV/c 2 and V = 0.05 GeV. 
The curves show the results of attempts to produce a 
peak in the vicinity of the data with interfering S-, P- 
and D-waves in the Kn~ channel. 

with the S-, P-, and D-wave intensities as shown in 
Fig. [5] Good fits to the data are obtained. The study 
of the angular dependence shows a clear asymmetry in 
cos 9k, where 6 k is the angle between the Kaon and 
ifj meson in the Kit~ rest frame (see Figs. 12 and 13 
in Ref. 0). 

The BaBar mass resolution at the Z(4430)~ mass 
is 7 (4) MeV/c 2 for the modes with J/ip (ip') (see 
Fig. 6 in Ref. [5]). Each event is efficiency corrected, 
and the conclusions are based on distributions ob- 
tained after background subtraction. The average effi- 
ciency varies slightly as a function of the Kit~ mass 
but can differ significantly between the different tp de- 
cay modes (see Fig. 7 in Ref. 0), however there are 
significant variations in efficiency with cos 6k for cer- 
tain regions of K-ir^ mass (see Fig. 36 in Ref. ll5l"). 

The comparison between the i/jit~ mass distribu- 
tion and the relevant Kir~ reflection was performed 
(see section 9 in Ref. [5]). This comparison is es- 
sential to investigate the need for a Z(4430)~ signal 
above the Kir~ reflection. Following the five Ktt~ 
regions (see section 10.D in Ref. (H) defined in the 




Figure 4: From BaBar Q, the m 2 ^. versus rn? K ^_ 
Dalitz plots for the final samples in the decay modes 

(a) B~ -»■ J/tpTr-K%, (b) B° -> J/t/m-K + , (c) 
B~ -> i//it-K° and (d) B° -»• i//ir~K + . 
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Figure 5: From BABAR [5], the fit results to the Ktt~ 
mass distributions for the combined Kit~ charge con- 
figurations (a) B~° ->• J/4m-K°> + and (b) B~'° -> 
i/j'tt~K 0,+ . The open dots show the data, and each 
curve represents different fit contribution, as indi- 
cated. 

Belle analysis (31, the projection is implemented in 
each Kir~ interval as shown in Fig. [6] In each fig- 
ure, a comparison (not a fit) between the data and the 
Ktt~ reflection is shown. Note that only one normal- 
ization factor for the Kn~ regions is used for each 
tp mode. The Kir~ reflections onto the ipn~ mass 
distributions describe well the data, with a small sta- 
tistical fluctuation at m K7T - < 0.795 GeV/c 2 . The 
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^7r~ mass distribution and the overall Kir~ projec- 
tion are shown in Fig. |7J No significant evidence 
for Z(4430)~ has been obtained. The Ktt~ projec- 
tion onto the tpir~ mass distribution is used as a non- 
peaking background. With such a background shape 
and a BW function for a Z(4430)~ signal, a signal is 
obtained for the overall Ktt~ region, but with shifted 
mass value. A similar small signal is obtained in the 
K*(892) and iff (1430) combined regions. A ~ 2a 
peak is obtained when the K*(892) and K|(1430) 
regions are vetoed. No (or a negative) signal for 
Z(4430)~ -> J/ipTT~ was found. 
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Figure 6: From BABAR [5], the ijm mass distri- 
butions in different regions of Kir~ mass for (a-e) 
B- fi -> J/ipTr~K ' + , and (f-j) £~ -> V'vr~K°' + ; 
the open dots represent the data, and the solid curves 
and shaded bands are due to the Kir~ reflections 
in the different i/j7r~ mass intervals, using the same 
overall normalization constant. In (f-j), the dot- 
dashed curves are obtained using Kn~ normalized 
moments for B~'° — > J/ipn~ K 0,+ , instead of those 
from B~'° — > ip'ir~ K°' + ; the vertical lines indicate 
m^r- = 4.433 GeV/c 2 . 

Although the Belle Collaboration used a data sam- 
ple equivalent to 1.46 times that of BABAR's, be- 
cause of different selection criteria BaBar has more 
events per fb _1 . A direct comparison between the 
two rriM W - distributions is shown in Fig. [8] For com- 
parison purposes, the BaBar distribution is not cor- 
rected for efficiency, however both Belle and BaBar 
distributions are background subtracted and shown in 
Fig. Uta) and Fig. Htb), respectively. The BaBar dis- 
tribution is scaled up by the factor 1.18 in order to 
normalize to the Belle data between the crosshatched 
regions of Fig. [He). The difference plot is shown in 




Figure 7: From BaBar [5], the i/jtt mass distribu- 
tions for the combined decay modes (a) B~'° — > 
J/ifjTT-K°'+ and (b) £~ -> ^'tt-K°' + . The points 
show the data (integrated over all Ktt~ regions) af- 
ter efficiency correction and background subtraction. 
The dashed curves show the Ktt~ reflection for a flat 
cos 9k distribution, while the solid curves show the 
result of cos 6k when accounting for the angular de- 
pendence (see Sec. 9 in Ref. Q). The shaded bands 
represent the effect of statistical uncertainty on the 
normalized moments. In (b), the dot-dashed curve 
indicates the effect of weighting with the normalized 
J/i/jtt~ K moments. The dashed vertical lines indi- 
cate the value of m^ n - = 4.433 GeV/c 2 . In (c) and 
(d), we show the residuals (data-solid curve) for (a) 
and (b), respectively. 



Fig. [He) where errors have been combined in quadra- 
ture; here x 2 /NDF = 54.7/58, which indicates that 
the two samples are statistically equivalent. 

Belle's Dalitz analysis of B — > t^'n~K 

After the BaBar Collaboration did not confirm the 
Z(4430)~ — > i/j'tt~ mass peak in their analysis of 
B — > Tp'-K^K decays (H, the Belle Collaboration per- 
formed a reanalysis [9] of their data that took detailed 
account of possible reflections from the Kit~ channel. 
Specifically, they modeled the B — > ip'ir^K process 
as the sum of two-body decays B — > tp'K* , where K* 
denotes all of the known K* -^Kir~ resonances that 
are kinematically accessible, and both with and with- 
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Figure 8: From BaBar Q, (a) The Belle dis- 
tribution after background subtraction; (b) the equiv- 
alent BaBar rriM n - distribution obtained after scal- 
ing due to different luminosities, (c) The difference 
between the mM w - distributions. The crosshatched 
area represents the exclusion region for the normal- 
ization procedure. The vertical-dashed line indicates 
m , v _ = 4430 MeV/c 2 . 



out a B — > Z(4430)~ K component. The results of 
this analysis, confirm the basic conclusions of Belle's 
initial publication (3J. 

Figure [9] shows the t/j'ir~ mass distribution in dif- 
ferent Kit~ intervals as obtained from the Dalitz-plot 
analysis by the Belle Collaboration [9]. In the region 
between the K*(892) and the # 2 *(1430) (0.982 < 
rriR-K < 1-332 GeV/c 2 , Fig. Etc)), a Z(4430)- peak 
is present. The data points in Fig. [10] show Af? /7r _ 
Dalitz plot projections with the prominent K* bands 
removed (as in Fig. [T]) compared with the results 
of the fit with no Z(4430)~ resonance, shown as a 
dashed histogram, and that with a Z(4430)~ reso- 
nance, shown as the solid histogram. The fit with the 
Z(4430)~ is favored over the fit with no Z(4430)" by 
6.4a. The fitted mass, M = 4443±J| t\l MeV/c 2 , 
agrees within the systematic errors with the earlier 
Belle result; the fitted width, T = 107±fjj Igg MeV, 
is larger, but also within the analysis' systematic er- 



rors of the previous result. In the default fit, the 
Z(4430)~ resonance was assumed to have zero spin. 
Variations of the fit that included a J = 1 assignment 
for the Z(4430)~ as well as models that included ad- 
ditional, hypothetical K* — ^AV - resonances with 
floating masses and widths, and radically different pa- 
rameterizations of the Ktt~ 5-wave amplitude do not 
change the conclusions. The product branching frac- 
tion from the Dalitz fit: B(B° -> Z(4430)~#) x 
£(Z(4430)~ V'vr") = (3.2j£t x 10" 5 is 
not in strong contradiction with the BaBar 95% C.L. 
upper limit of 3.1 x 10" 5 0. 
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Figure 9: From Belle I0, M?,^_ projections for Mktv 
bands: (a) below the K* (892); (b) at the if* (892); 
(c) between the #*(892) and the #£(1430); (d) at 
the #|(1430); and (e) above the # 2 *(1430). The dots 
show the data, the solid (dotted) histograms are the re- 
sults for the model with (without) a single tp'ir~ state, 
and the dashed histograms are the background con- 
taminations. 
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Figure 10: From Belle (9J, the M\,^_ projection of 
the Dalitz plot with the K* bands removed is shown as 
data points. The histograms show the corresponding 
projections of the fits with and without a Z(4430) _ — >■ 
ip'ir~ resonance term. 



Belle's two Z peaks in the Xa^ channel 

In addition to the Z(4430)~, Belle has presented 
results of an analysis of B — > Xci^ ~K decays 
that require two resonant states in the Xci^r chan- 
nel 0]. The M 2 _ vs. M 2 _ Dalitz plot, shown 
in Fig.[TT] shows vertical bands of events correspond- 
ing to K*(892)->Kir- and K* ( 1430) ^Kir~, plus a 
broad horizontal band near M 2 _ ~17.5GeV 2 /c , 
indicating a possible resonance in the Xci^r chan- 
nel. In this case, this horizontal band corresponds to 
cos 8 n ~ ~ 0, a location where interference between 
partial waves in the Kit~ channel can produce a peak 
and, thus, a detailed Dalitz analysis is essential. 

For B — > Xcin~K, the kinematically allowed 
mass range for the Ktt~ system extends beyond the 
1^3(1780) F-wave resonance and S-, P-, D- and 
F-wave terms for the Kn~ system are included in 
the model. The fit with a single resonance in the 
Z~ — > Xci^r channel is favored over a fit with only 
K* resonances and no Z~ by more than 10<r. More- 
over, a fit with two resonances in the Xci^~ channel 
is favored over the fit with only one Z~ resonance by 



Figure 11: From Belle [41, the M\ _ (horizontal) vs. 
M 2 ivr _ (vertical) Dalitz plot distribution for candi- 
date B — > Xci^K events. 



5.7a. The fitted masses and widths of these two reso- 
nances are: Mi = 4051 ± 14^° MeV/c 2 and r x = 



,+21 +47 



-17 



% MeV and M 2 = 4248 



+44 +180 
-29 -35 



MeV/c 2 



177^o t 3 6 \ 6 MeV. Here also, variations 



and T 2 - x. ,_ 39 
of the fit that use a J = 1 assignment for the Z~ 
states and models that include additional, hypotheti- 
cal K* -^Kir~ resonances with floating masses and 
widths, and different parameterizations of the Kir~ 
5-wave amplitude do not change the conclusions. 

The product branching fractions for B — > Z~(—> 
Xci^ )K have central values similar to that for the 
Z(4430)~ but with large errors. Figure (T2] shows 
the M Xcl7r - projection of the Dalitz plot with the 
K* bands excluded and the results of the fit with no 
Z~ — > Xci^r resonances and with two Z~ — > Xci^r 
resonances. 

Summary 

If the peaks reported by Belle in the t^'ir" and 
Xci7T~ channels are in fact meson resonances, they 
would be "smoking guns" for exotics. It is there- 
fore important that the Belle results are confirmed 
(or refuted) by other experiments. Both the BABAR 
and Belle ip'n~ analyses could benefit from additional 
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Figure 12: From Belle [4], the data points show the 
M Xcl7T - projection of the Dalitz plot with the K* 
bands removed. The histograms show the correspond- 
ing projections of the fits with and without the two 
Z — > Xcitt resonance terms. 

statistics. However, both experiments have completed 
data taking and the BaBar analysis used their com- 
plete dataset while the Belle analysis used about 75% 
of their final total data sample. Approximately 3 ab _1 
of integrated luminosity would be necessary to bring 
the ip'7r~ sample to the statistical level of BaB^R's cur- 
rent J/tpTT" sample, and this amount of data will not 
be available until the operation of one of the Super- 
B factories, which will not happen for at least a few 
years. In the meantime, it would be valuable if BaBar 
performed a search for the Z\ (4050) ~ and Z 2 (4250) ~ 
in the Xci^r channel and Belle reproduced BaBar's 
study of the J/i/jtt~ channel using the existing data 
sets. The DO and CDF experiments at the Tevatron 
have data on hand that could be used to carry out 
searches for inclusive production of the Z(4430)~ 
and we hope that results will be available from them 
in the near future. 
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